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ABSTRACT 

We describe the analysis of the seven broad-band X-ray continuum observations of the archetypal Seyfert 1 galaxy 
NGC 5548 that were obtained with XMM-Newton or Chandra, simultaneously with high-energy (> 10 keV) observations 
with NuSTAR and INTEGRAL. These data were obtained as part of a multiwavelength campaign undertaken from the 
summer of 2013 till early 2014. We find evidence of a high-energy cut-off in at least one observation, which we attribute 
to thermal Comptonization, and a constant reflected component that is likely due to neutral material at least a few 
light months away from the continuum source. We confirm the presence of strong, partial covering X-ray absorption 
as the explanation for the sharp decrease in flux through the soft X-ray band. The obscurers appear to be variable in 
column density and covering fraction on time scales as short as weeks. A fit of the average spectrum over the range 
0.3 - 400 keV with a realistic Comptonization model indicates the presence of a hot corona with a temperature of 401 jq 
keV and an optical depth of 2.71° 2 ^ spherical geometry is assumed. 

Key words, galaxies: active — galaxies: Seyfert — X-rays: galaxies 


1. Introduction 


The high-energy spectrum of active galactic nuclei (AGNs) 
is typically dominated by a primary power law component, 
which is generally believed to originate in a hot plasma, the 
so-cal led corona, localized in the inner part of the accretion 
flow (Reis & Millerf 20131. The presence of a high-energy 
cut-orf around ~100 keV has been observed in several 


sources (see, e.g., |Perola et al. 

2002 

Malizia et al. 

2014 

Brenneman et aL||2014| |Marinucci et aL||2014| |Ballantyne 

et al.|2014) and is commonly int 

srpreted as the signature of 
y, the optical/UV photons 

thermal Comptonization. Name 


from the accretion disk are upscattered by the ho t electrons 
of the corona up to X-ray energies (see, e.g., Haardt & 


Maraschi 1991 


emission can be modibecl by different processes 
absorption from surrounding gas and Compton reflection 
from the disk itself or more distant material. Moreover, the 
spectra of many AGNs rise smoothly below 1-2 keV above 


Haardt et 
modified 


al. 

bF 


1994 19971. This 


different 


primary 
such as 


the extrapo lated high-energy power law (see, e.g., Bianchi 
et al. 1(2009 ). The nature of this so-called soft excess is still 


under debate (see, e.g.. Done et al. ][m^ . 


The geometrical and physical properties of the differ¬ 
ent structures leading to the different spectral components 
are still poorly known. First and foremost, a more com¬ 
plete understanding of the high-energy emission of AGNs 
requires disentangling the spectral properties of each com- 
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ponent and constraining their characteristic parameters. 
For this goal, the availability of multiple broad-band 
observations with a high signal-to-noise ratio is essential. 
Using the data from a previous long multiwavelength 
campa i gn on Mrk 509 (Kaastra et al. 2011), |Petrucci 
et al. (2013) made a detailed study of the high-energy 
spectrum of that AGN. In that case, the simultaneous 
data covering the range from optical/UV to gamma-ray 
energies made it possible to test Comptonization models 
in full detail and provide a physically motivated picture of 
the inner regions of the source. The hard X-ray continuum 
was well fitted by a Comptonized spectrum produced 
by a hot, optically thin corona with kT^. ~ 100 keV and 
T ~ 0.5, proba bly localized in the ve ry inner part of the 
accret ion flow (Petrucci et al. 2013). Mehdipour et al. 
(2011), studying the optical to A-ray variability, found a 
correlation between the optical/UV and soft (<0.5 keV) 
X-ray flux, but no correlation between the optical/UV and 
hard (>3 keV) X-ray flux. This result was interpreted as 
indicating the presence of a second, "warm" corona with 
kTp ~ 1 keV and r ~ 15, responsible for the soft excess as 


well as the optical/ UV emission (Mehdipour et al. 2011 
Petrucci et al.|[2M3 l. 


In the present paper, we discuss results based on a multi¬ 
wavelength monitoring program on the AGN NGG 5548, 
an X-ray bright Seyfert 1 galaxy, h osting a supermassive 


and h 

y^ing at 

et al. 

1991 


as given in the NASA/IPAG Extragalactic 


.in¬ 


ability typica]_ofan_u^bscured AGN of that black hole 


mass (Ponti et al. 2012). High spectral resolution X-ray 
and UV observations of this source have previously shown 
a typical ionized outflow, seen as the "warm absorber" 
through narrow lines absorbing the UV and soft X-ray 
emission from the core (|Crenshaw et al.||2003l |Steenbrugge 


The multiwavelength campaign was conducted from 
May 2013 to February 2014, involving observations with 
the XMM-Newton, HST, Swift, NuS TAR, INTEGRAL and 
Chandra satellites (Kaastra et al. 2014 hereafter K14). 
During the campaign, the nucleus was found to be obscured 
by a long-lasting and clumpy stream of ionized gas never 
seen before in this source. This heavy absorption manifests 
itself as a 90% decrease of the soft X-ray continuum 
emission, associated with deep and broad UV absorption 
troughs. The obscuring material has a distance of only a 
few light days from the nucleus, and it is outflowing with 
a velocity up to five time s faster than those observed in 
the warm absorber (K14l. Given its properties, this new 
obscuring stream likely arises in the accretion disk. 

This work is part of a series of papers focusing on different 
aspects of the campaign. Here we concentrate on the 
broad band (0.3-400 keV) X-ray emission, focusing on 
the observations featuring high-energy (>10 keV) data 
provided by NuSTAR and INTEGRAL. 

The paper is organized as follows. In Sect. we de¬ 
scribe the observations and data reduction. In Sect. |3] 
we present the analysis of the hard (>5 keV) X-ray and 
total (0.3-400 keV) X-ray spectrum, fitted with both a 
phenomenological and a realistic Gomptonization model. 
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Fig. 1. Light curves in different energy ranges. On the x-axis, we 
report the last three digits of the Modified Julian Date. On the y- 
axis, on both left and right side, we report the count-rate. Black 
circles: XMM-Newton (0.3-0.5 keV, 0.5-2 keV, 2-10 keV); ma¬ 
genta nablas: Chandra (0.5-2 keV, 2-10 keV); red squares: NuS¬ 
TAR (20-40 keV, 40-80 keV); blue triangles: INTEGRAL (20-40 
keV, 40-80 keV, 80-300 keV). The vertical dotted lines highlight 
the seven observations we analysed in this paper (i.e., those in¬ 
cluding NuSTAR and/or INTEGRAL data). 


In Sect. we discuss the results and summarize our 
conclusions. 


2. Observations Data Reduction 

The multi-satellite campaign was led by XMM-Newton, 
which observed NGG 5548 twelve times during 40 days 
in June and July 2013, plus once in December and 
once in February 2014. For the full description of the 
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Table 1. The logs of the simultaneous XMM-Newton, NuSTAR and/or INTEGRAL observations of NGG 5548 during our cam¬ 
paign. 


Obs. 

Satellites 

Obs. Id. 

Start time (UTG) 
yyyy-mm-dd 

Net exp. 
(ks) 

1 

XMM-Newton 

0720110401 

2013-06-30 

38 


INTEGRAL 

10700010001 

2013-06-29 

62 

2 

XMM-Newton 

0720110601 

2013-07-11 

37 


NuSTAR 

60002044002/3 

2013-07-11 

50 


INTEGRAL 

10700010002 

2013-07-11 

50 

3 

XMM-Newton 

0720110701 

2013-07-15 

37 


INTEGRAL 

10700010003 

2013-07-15 

50 

4 

XMM-Newton 

0720111101 

2013-07-23 

38 


NuSTAR 

60002044005 

2013-07-23 

50 


INTEGRAL 

10700010004 

2013-07-23 

52 

5 

Chandra 

16314 

2013-09-10 

120 


NuSTAR 

60002044006 

2013-09-10 

50 

6 

XMM-Newton 

0720111501 

2013-12-20 

38 


NuSTAR 

60002044008 

2013-12-20 

50 

7 

XMM-Newton 

0720111601 

2014-02-04 

38 


INTEGRAL 

11200110001 

2014-01-17 

94 



11200110002 

2014-01-22 

40 



11200110003 

2014-02-09 

30 


program, see K14 and Mehdi pour et al. ([submit ted I. 
During the ca mpaign, NuBTAH ([Harrison et al.||2013 1 and 
INTEGRAL (Winkler et al. 2003|) regularly observed the 
source simultaneously with XMM-Newton. A simultaneous 
Chandra/NuSTAR observation was also performed in 
September 2013. The log of the data sets analysed in 
this paper is reported in Table with their observation 
ID, their start date, and net exposure time. The latest 
XMM-Newton observation is dated February 4, and the 
corresponding spectrum has been analysed jointly with an 
average spectrum obtained from three INTEGRAL expo¬ 
sures obtained in January and February, which show no 
significant variations. 


tribution as of version 6.14), using calibration files from 
NuSTAR CALDB v20130710. The raw event files were 
cleaned with the standard depth correction, which greatly 
reduces the internal background at high energies, and 
data from the passages of the satellite through the South 
Atlantic Anomaly were removed. Spectra and light curves 
were extracted from the cleaned event files using the 
standard tool NUPRODUCTS for each of the two hard X-ray 
telescopes aboard NuSTAR, focal plane modules A and 
B (FPMA and FPMB); the spectra from these telescopes 
are analysed jointly in this work, but are not combined. 
Finally, the spectra were grouped such that each spectral 
bin contains at least 50 counts. 


In all the XMM-Newton observations, the EPIC in¬ 
struments were operating in the Small-Window mode 
with the thin-filter applied. For our analysis, we used the 
EPIC-pn data only. The data were processed using the 
XMM-Newton Science Analysis System (SAS vl3). The 
spectra were grouped such that each spectral bin contains 
at least 20 counts. Starting in 2013, a gain problem with 
the pn data appeared, causing a shift of about 50 eV near 
the Fe-K complex (see the detailed discussion in Cappi et 
al., in prep.). In the present paper this shift was corrected 
for by leaving the redshift free to vary in our models. 
Moreover, this small offset yields poor fits near the energy 
of the gold M-edge of the telescope mirror, and for this 
reason the interval between 2.0 and 2.5 keV was omitted 
from the fit. However, this problem has no effect upon the 
broad-band modelling, in particular the Fe-K line strength 
is correct, thus allowing for the reflection component to be 
correctly modelled. 

The NuSTAR data were reduced with the standard 
pipeline (nupipeline) in the NuSTAR Data Analysis 
Software (nustardas, vl.3.1; part of the heasoft dis- 


INTEGRAL data were reduced with the Off-line Sci- 


entiflc Analysis software OSA 10.0 (Courvoisier et al. 


20031. IBIS/ISGRI spectra have been extracted using the 


standard spectral extraction tool ibis_science_analysis 
and combined together using the OSA spe_pick tool 
to generate the total ISGRI spectrum, as well as the 
corresponding response and ancillary files. 


The Chandra HRG/LETG spectra (first order) of the 
120 ks observation of September 10 were reduced with the 
package CIAO v4.6 and reprocessed with the standard tool 
chcuidr a_repr o. 


Finally, we made use of archival BeppoSAX and Suzaku 
data to compare the results of the campaign with past 
X-ray observations of the source. We focused on a 1997 
BeppoSAX long (^300 ks) observation and on seven 
Suzaku observations (~30 ks each) performed in 2007. 
For a description of the observations and data reduction 
procedure we refer tojNicastro et al.|(2000l for BeppoSAX, 
and Liu et al. (20101 for Suzaku. 
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Fig. 2. Upper panel: the 10-79 keV spectra of the four NuS- 
TAR observations. Lower panel: residuals when simultaneously 
fitting with a power law. For clarity, we used only FPMA data, 
which are rebinned in the plot. 


In Fig. [T] we show the light curves of all the four¬ 
teen XMM-Newton observations, the four NuSTAR 
observations and the five INTEGRAL observations of the 
campaign. The variability shown by the light curves in the 
low-energy ranges is mainly due to obscuration, as shown 
in Sec. We will also show that the primary emission 
exhibits intrinsic variations as well, which explain the 
va riability observed a t higher ene rgies. We refer the reader 
to Mehdipour et al. (submittedI for the full time-line of 
the campaign (see their Fig. 1). 


3. Spectral analysis 


Spectral analysis and model fitting w as carried out with the 
XSPEC 12.8 package (Arnaud 19961, using the minimi¬ 
sation technique throughout. In this paper, all errors are 
quoted at the 90% confidence level. 


3.1. Spectral curvature above 10 keV 

In Fig. we plot for comparison the four NuSTAR/FI’M.A 
spectra above 10 keV, simultaneously fitted with a power 
law with tied parameters (x^/dof = 1549/1178). The pho¬ 
ton index obtained is F = 1.58 ± 0.02. The plot shows that 
the spectral variability in the hard X-rays during the cam¬ 
paign isnot dramatic, while it is prominent in the soft X- 
rays (K14 Cappi et al., in prep.; di Gesu et al., in prep.) 
where obscuration plays a major role. After getting a gen¬ 
eral idea about the power law parameters, we tested the 
presence of spectral curvature by individually fitting each 
of the four NuSTAR spectra above 10 keV, first with a sim¬ 
ple power law and as a second step including an exponential 
cut-off. The addition of the cut-off results in an improve¬ 
ment in all fits (Ay^ ~30), with a probability of chance 
improvement less than 10“® according to the F'-test. This 
result motivates the use of a cut-off power law in the analy¬ 
sis of the data. However, as a next step we test whether the 
spectral curvature is partly due to a reflection component. 
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3.2. The reflection component 


To test for the presence of a reflection component, we fitted 
the data above 5 keV. By choosing 5 keV as a lower energy 
limit, we avoid the effects due to the strong absorption 
in the soft band. Good statistics are needed, so we only 
consider the four observations including NuSTAR data 
(obs. 2, 4, 5 and 6 of Table IB. 

We fitted the XMM-Newton/pn, NuSTAR and INTE¬ 
GRAL data simultaneously, leaving the cross-calibration 
constants free to vary. The normalization factors with 
respect to pn were always between 1 and 1.1. In partic¬ 
ular, the spectral agreement between XMM-Newton and 
NuSTAR is generally good in their common 3-10 keV 


bandpass (Walton et al. 2013 20141. We used a simple 


consisting of a cut-off power law modified by 
plus reflection. We included a neutral ab- 


model, 
absorption, 

sorber with covering fraction fixed to 1, while the column 
density was free to vary. This is clearly an arbitrary 
choice, but the way we model the absorption, with a lower 
energy limit of 5 keV, does not affect the results that follow. 


Following the analysis done by Gappi et al. (in prep.) 
and motivated by their detection of a neutral, narrow 
and constant Fe Ka line, we used the PEXMON model 
(Nandra et al. 20071 in XSPEC for the reflection compo¬ 
nent. PEXMON combines a neutral Gompton reflection, as 


described by the pexrav model (Magdziarz & Zdziarski 


19951 in XSPEC, with Fe Ka, Fe K/t, JNi Ka lines and the 


F'e Ka Gompton should er, generated self-consi stently by 
Monte-Garlo simulations (George & Fabian|1991 ). We fixed 
the inclination angle of the reflector to 30 d eg, appropriate 
for a type 1 source (e.g., Nandra et al. ]I9^ and consistent 


with past analyses of the Fe K lines in NGC 5548 (Yaqoob 


et al. 


2001 Liu et al. 20101, but also with the inclination 


estimated for the broad-line region (Pancoast et al. 


and for the disk-corona system (Petrucci et al. 20001. The 


20141 


photon index, cut-off energy and normalization of the 
incident radiation were free to vary and not tied to those of 
the primary power law. This choice is justified a priori by 
the observed constancy of the Fe Ka line, and a posteriori 
by our constraints on the reflection component, which is 
consistent with being constant during the campaign. In this 
case, the distance of the reflecting material from the X-ray 
emitting region should be at least a few light months, i.e. 
a few X 10“^ gravitational radii, as we will discuss in Sect.|^ 


The 90% contour plots of cut-off energy and normal¬ 
ization vs. photon index of PEXMON are plotted in Figs. 
and 1^ The photon index and cut-off energy of the reflected 
continuum are not well constrained, however the contours 
are mutually consistent. The normalization vs. photon 
index contours are very similar for the four observations, 
with indication of a slightly larger normalization for obs. 
6; nevertheless they are also marginally consistent with 
each other. These results are consistent with the reflected 
emission coming from remote material, and thus the 
parameters being constant during the campaign. For this 
reason, we chose to freeze hereafter the parameters of the 
reflection component to their weighted average derived 
from our analysis. We thus fixed the PEXMON photon 
index to 1.9, cut-off energy to 300 keV and normalization 
to 5.5 X 10“^ (xsPEC units). The photon index of 1.9 is 
relatively high compared to the average value for NGG 
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Fig. 3. PEXMON cut-off energy vs. photon index 90% confidence 
level contours for the 4 NuSTAR observations, fitted above 5 
keV. The grey solid lines indicate the values we chose to fix in 
the later broad-band analysis (1.9 for the photon index, 300 keV 
for the cut-off energy). 



1.4 1.6 1.8 2 2.2 


PEXMON Photon Index 

Fig. 4. PEXMON normalization vs. photon index 90% confidence 
level contours for the 4 NuSTAR observations, fitted above 5 
keV. The grey solid lines indicate the values we chose to fix in the 
later broad-band analysis (1.9 for the photon index, 5.5 x 10“® 
for the normalization). 


3.3. Broad-band spectral analysis 

Having obtained constraints on the parameters of the re¬ 
flection component using the data above 5 keV, we now 
proceed to fit the total X-ray energy range available. We 
first use a phenomenological model, which is described be¬ 
low. Given the different satellites used in the seven obser¬ 
vations, the energy range is limited to 0.5-79 keV for obs. 
5 {Chandra/NuSTAR) and 0.3-79 keV for obs. 6 {XMM- 
Newton/NuSTAR). For all the other observations, the en¬ 
ergy range is 0.3-400 keV, thanks to the addition of the 
INTEGRAL data. 


3.3.1. The model components 

Following Sect. |3.1| and |3.2[ we modelled the continuum 
with a cut-off power law and the reflected component with 
PEXMON, whose parameters were fixed. We describe in the 
following the other spectral components. 


Obscuration: The model presented by |K14 includes two 
ionized, partially covering absorbers. Both were found to 
have a low degree of ionization, namely log^ < —4 and 
log^ ~ —1-2, such that at the EPIC-pn energy resolution 
they are almost equivalent to neutral absorbers. In partic¬ 
ular, the absorber with log^ = —1. 2 was found to have a 
column density of 1.2 x 10^^ cm“^ (K14l and such an ab¬ 
sorber is indistinguishable, using pn data only, from a neu¬ 
tral one with a column density of 1.1 x 10^^ cm“^. Since a 
detailed modelling of the obscuration is out of the scope of 
the present paper, we modelled the two absorbing compo¬ 
nents with PCFABS in XSPEC, which describes partial cov¬ 
ering absorption by a neutral medium. Both the column 
density and the covering fraction of the two components 
were free to vary between different observations. 


The warm absorber: The model by K14| includes six 


warm absorber components, with different column densities 
{Nfj ~ 10^°-10^^ cm“^), ionization parameters (log^ ~ 1- 
3), and outflow velocities {v ~ 250-1000 km/s). We pro¬ 
duced a table, with fixed parameters, from the best-fit 
model of the warm absorber derived by |K14| and included 
it in XSPEC. For a comprehensiv e analysis o f the UV pr op- 
erties of the warm absorber, see Arav et al. (submitted I. 


5548, which is ~ 1.7 (see, e.g., Sobolewska fc Papadakis 
20091 but has indeed been observed in the pasFlNandra 


et al 


1991 Reynolds 1997, Chiang et al. 2000). More 


importantly, the Swift monitoring discussed in Mehdipour 
et al. (in prep.) shows that the photon index reached 
1.8-1.9 during most of 2012. Finally, from the PEXMON 
parameters we estimate the average reflection fraction / as 
0.3 ±0.1, calculated as the ratio between the 20-40 keV flux 
illuminating the reflecting material and the 20-40 keV flux 
of the primary cut-off power laT] For a distant reflector, 
this ratio might not have a clear link with the covering 


facto r as seen from the X-ray source ( Malzac fc_Petrucci 
2002), but this is discussed in more detail in Sec.|§ 


^ This definition is slightly different from the traditional R — 
n/27r, where Q is the solid angle subtended by the reflector 
(Magdziarz fc Zdziarski|1995 1. 


The soft X-ray excess: The physical origin of this compo¬ 
nent is still under debate, and it can be modelled with sev¬ 
eral different spectral models. The most common interpre- 
tations rely on relativistically blurred reflection (see, e.g.. 


Crummy et al. 2006 Cerruti et al. ]ml| [Walton et al.|2012t 


or Comptonization by a warm mediuni (see, e.^ Petrucci 

et al. 12013 [Boissay et al.|[2M4 Matt et Sl|2014 ). Detailed 
modelling of the soft X-ray excess is beyond the scope of this 
paper, we thus followed the analysis done by Mehdipour et 
al. (in prep.) and assumed a thermal Comptonization com¬ 


ponent modelled by COMPTT (Titarchuk 
The parameters of COMPTT were frozen, 


1994 ) in XSPEC. 
for each individ¬ 


ual observation, to the values reported by Mehdipour et al. 
(in prep.), to which we refer the reader for details. We note 
that Mehdipour et al. (in prep.) found a significant flux 
variability of the soft excess between different observations, 
which we take into account here, but no apparent spectral 
shape variability. 


Article number, page 5 of 


10 









































































A&A proofs: manuscript no. articolo 




Photon Index 


Fig. 5. Data and best-fit model for obs. 2 (XMM- 
Newton/NuSTAR/INTEGRAL, see Table [^. Upper panel: 
broad-band data and folded model. Middle panel: contribution 
to ■ Lower panel: best-fit model f{E) with the plot of the 
various additive components, namely the cut-off power law, the 
reflection component, the soft excess and the narrow lines. 


Narrow emission lines: A number of emission lines, such as 
[O VII] at 0.56 keV and [C V] at 0.3 keV, were detected in 
the XMM-Newton RGS spectra, as shown by |K14| Again, 
we produced a table from the best-fit model and included it 
in XSPEC. These lines arguably originate from the Narrow- 
Line Region, therefore we do not expect lin e flux variations 
during the campaign (Peterson et al.pOlSl. 


3.3.2. Best-fit results 

The best-fit parameters are reported in Table while 
in Fig. 1^ we show the data, residuals and best-fit model 
for obs. 2 {XMM-Newton/NuSTAR/INTEGRAL). The 
broad-band fit allows us to constrain the spectral shape 
more reliably with respect to the fit above 5 keV. In Fig. 

and we show the contour plots of the high-energy 
cut-off and normalization versus photon index, both for 
the primary continuum. The contours show significant 
variability of the photon index from ^ 1.5 (obs. 6, 7) to 
^ 1.7 (obs. 1, 2, 4, 5). It should be noted that the photon 
index given in Kaastra et al. (2014), i.e. F = 1.566 ± 0.009, 
is a mean value which comes from the analysis of the 
co-added XMM-Newton and NuSTAR spectra of summer 
2013. By including the individual NuSTAR and INTE¬ 
GRAL high-energy spectra, we could constrain the photon 
index for each observation. Furthermore, the cut-off energy 
shows significant variability as well. We derive a lower limit 
for the cut-off energy in all the observations, and a measure 
of 70l^° keV for obs. 6. This measure is inconsistent, at 
90% confidence level, with the lower limits determined 
from the other observations, with the exception of obs. 3. 
Finally, we note that the constraints we derived on the 
cut-off energy are mostly compatible with the average 
value estimated in sec. |3.2| for the reflection component 
alone, i.e. 300 keV. 

The broad-band fit helps constrain the obscurer com¬ 
ponents as well. In Fig. and we plot the contours of 
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Fig. 6. Cut-off energy vs. primary photon index 90% confidence 
level contours for obs. 1 (black), obs. 2 (red), obs. 3 (green), obs. 
4 (blue), obs. 5 (cyan), obs. 6 (magenta), obs. 7 (orange). 


the covering fraction versus the column density of the two 
neutral absorbers, while in Fig. [T^ we plot the contours 
of the two column densities versus each other. These 
contours show that absorption variability is statistically 
significant for both the low and high-Np absorbers. 
This can be seen, for example, by comparing obs. 2 (red 
contours) with obs. 4 (blue contours). Since the time 
interval between obs. 2 and 4 is ~ 10 days, we can conclude 
that the variability is significant on this time scale. In two 
cases, namely obs. 5 and 7, only one, low-Ab/ absorber 
is needed, and adding the high-Np component does not 
improve the fit. For obs. 5 this can be explained by the 
lower sensitivity of Ghandra HRC/LETG with respect to 
XMM-Newton EPIG-pn, which results in a lower S/N in 
the soft band. For obs. 7, which is the last one and dated 
February 4, the absence of the high-Np component could 
simply mean that the source started changing back to 
the more usual, unobscured condition. However, for the 
model described above, we only obtain a relatively poor fit 
(x^/dof = 651/424), with strong, positive residuals around 
0.5 keV, and negative residuals around 0.8 keV. We stress 
that the models we used for the warm absorber and the 
narrow emission li nes w ere obtained from the analysis of 
pn and RGS data (K14), while our analysis is based on pn 
data alone. We may thus attribute such residuals in obs. 7 
to a combination of the previously mentioned gain problem 
of the pn data with cros s-calibration issues between the pn 
and RGS instruments (Kirsch et al. 20041. We obtained 
a much better fit (Ax^ ^ —200) by adding a narrow 
Gaussian line with GAUSS in XSPEC, and an absorption 
edge with edge, without affecting the parameters of the 
primary cut-off power law. The emission line was found 
at E = 0.54 ± 0.01 keV, with an upper limit on the 
equivalent width of 20 eV. The absorption edge was found 
at E = 0.80 ± 0.01 keV with r = 0.28 ± 0.06. 


3.4. Test of Comptonization model 

The last step of our analysis consists of fitting the average 
spectrum with a more realistic model, where the primary 
continuum is produced self-consistently via Gomptoniza- 
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Table 2. Best-fit parameters of the primary cut-off power law and the two absorbers (1) and (2), htting the data down to 0.3 keV. 
In the case of obs. 7, the model included a Gaussian line and an absorption edge (see text). 

Ec is in keV, Fi kev is in units of 10“^ ph/keV cm^ s and Nh is in units of 10^^ cm”^. 

I The second absorber was not needed in obs. 5 {Chandra plus NuSTAR) and obs. 7. In these two cases, fixing N^'’ to 8.5 resulted 
in = 0. 


Obs. 

r 

Ec 

±lkeV 

ivW 

7^(1) 


^(2) 

X^/dof 

1 

1 70+0.04 

‘0'_o.05 

>120 

0.96 ±0.09 

0.75 ±0.03 

0 7Q7+0-006 

u.<y/_o.oo7 

9± 1 

n 

'o-^^-o.oe 

538/424 

2 

1 Y^+O.Ol 
-L- <"±_0.04 

>250 

1 1 O+0.03 

-L-- 1 - 0 _o.o 7 

0.85 ±0.02 

O.829l°;004 

n 2+0.S 

n 09+0.02 

1587/1393 

3 


>50 

0.7+oa 

-L-UD_o.io 

O-S+Hs 

9^3 

n 97+0.06 
-0.10 

508/420 

4 

1.68 ±0.05 

>120 

o-87;n6 

1.27 ±0.06 

0.77 ±0.01 

8.5t°;8 

0.46 ±0.03 

1346/1329 

5 

1 74+0.01 
'^-0.03 

>180 

1.00 ±0.05 

1.8 ±0.02 

0.62 ±0.03 

t 

t 

1113/1094 

6 

1.49 ±0.05 

70+40 

<U_io 

0.56 ±0.05 

4 -^- 0 .1 

0.68 ±0.04 

8±2 

0.33 ±0.06 

1484/1308 

7 

1.49 ±0.02 

>120 

0.50 ±0.01 

U.++o.o4 

0.678 ±0.007 

t 

t 

452/420 


Table 3. Best-fit parameters of the average data set, using the realistic Comptonization model compps for the continuum. The 
physical units are the same as in Table 

The normalization is dehned as A = Rkm/Dio kpc where Rkm is the source radius in km and Dio kpc the distance in units of 10 
kpc. 


kTc (keV) 

y A(xlO^) 

ivW 


nP 

c‘» 

X^/dof 

0 0 

+ 1 

0 

0.85^°:'° 6.0;°;4 

0.99 ±0.02 

0.766 ±0.003 

9.6 ±0.5 

0.36 ±0.01 

3079/2627 



Photon Index 


Column Density (10^^ cm-^) 


Fig. 7. Normalization vs. primary photon index 90% conhdence 
level contours for obs. 1 (black), obs. 2 (red), obs. 3 (green), obs. 
4 (blue), obs. 5 (cyan), obs. 6 (magenta), obs. 7 (orange). 


Fig. 8. Obscurer component 1: Nh — Cf 90% confidence level 
contours for obs. 1 (black), obs. 2 (red), obs. 3 (green), obs. 4 
(blue), obs. 5 (cyan), obs. 6 (magenta), obs. 7 (orange). 


tion of soft photons emitted by the disk. A more detailed 
analysis of the entire campaign with realistic Comptoniza¬ 
tion models is deferred to a future work. 

First, we produced an average spectrum merging 
respectively the XMM-Newton, NuSTAR and INTE¬ 
GRAL spectra used in this paper (Table with the 
MATHPHA tool of HEASOFT. We then proceeded to fit 
this average spectrum over the range 0.3-400 keV, at first 
with t he phe nomenological model we have discussed above 
(Sect. 3.3.11. The resulting fit had y^/dof = 3314/2631, 


with strong, positive residuals around 0.5 keV, and 
negative residuals around 0.6 keV (analogously to the 
already discussed obs. 7). We thus included two narrow 
Gaussian lines, with E — 0.53 ± 0.06 keV (in emission, 
equivalent width < 50 eV) and E = 0.60 ± 0.02 keV (in 
absorption, equivalent width < 20 eV). Again, we note 
that such residuals are most likely un physical and due t o 
cross-calibration between pn and RGS ( Kirsch et al.]2004 ). 
However, the addition of these lines/ while improving 
the fit, does not affect the best-fit parameters of the 
cut-off power law, namely T = 1.67 ± 0.02 and Ec > 250 
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Fig. 9. Obscurer component 2: Nh — Cf 90% confidence level 
contours for obs. 1 (black), obs. 2 (red), obs. 3 (green), obs. 4 
(blue), obs. 6 (magenta). This obscurer was not found in obs. 5 
and 7. 


Fig. 11. kTs — y 90% confidence level contours for the spectra 
fitted with COMPPS, using spherical geometry. Black: average 
spectrum of the 2013 campaign; red: average spectrum of the 
seven Suzaku observations taken in 2007; green: BeppoSAK ob¬ 
servation taken in 1997. 



Fig. 10. Nh — Nh 90% confidence level contours for obs. 1 
(black), obs. 2 (red), obs. 3 (green), obs. 4 (blue), obs. 6 (ma¬ 
genta) . 


keV, with x^/dof=3080/2627 and no prominent residual 
features. The relatively large is most likely an effect 
of averaging variable spectra (in particular in the soft 
band, due to absorption variability), and of the not perfect 
cross-calibration between XMM-Newton and NuSTAR, 
which especially appears in the co-added spectrum. We 
also note that XMM-Newton and NuSTAR do not cover 
strictly simultaneous times, since the exposure time of the 
two satellites is slightly different (see Table . 


We then replaced the cut-o ff power law with the Comp- 
tonization model COMPPS (Poutanen & Svensson 19961 
in XSPEC. COMPPS models the thermal Comptonization 
emission of a plasma cooled by soft photons with a disk 
black-body distribution. We chose a spherical geometry 
(parameter GE 0 M =-4 in COMPPS) for the hot plasma and 
a temperature at the inner disk radius of 10 eV, fitting 
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for the electron temperature kT^ and Compton parameter 
y = AT{kTg/m(,c'^). The best-fit parameters are reported in 


(black contour). From the best-fit temperature kTe = 40 
keV and parameter y — 0.85, we derive an optical depth 
T = 2.7l5'2, where we estimate the errors from the kTg — y 
contour. 


Table ^ and the kT^, — y contour plot is shown in Fig. 11 


4. Discussion and conclusions 

The broad-band coverage of the XMM-Newton, NuSTAR 
and INTEGRAL data allowed us to study the different 
components of the high-energy spectrum of NGC 5548 in 
detail. 


First, we have been able to constrain the reflected 
component, which is consistent with being constant on 
a few months time scale, while the primary emission 
exhibits significant variability on a few days time scale. 
For example, the flux in the 2-10 keV band exhibits 
variations up to a factor of two over about 5 days during 
the campaign, with an outburst of several days du ration in 
August-September (Mehdipour et al. submitted see their 
Fig. 2). The lack of variability of the reflection component 
indicates that reflection is produced by cold and distant 
matter, lying at least a few light months away from the 
central X-ray source. This is consistent with the analysis 
of the Fe Ka line by Cappi et al. (in prep.) and with 
the absence of a significant cont ribution from relativisti c 
reflection in this source (see, e.g., Brenneman et al.||2012|. 
The lower limit of the reflector distance is consistenit with 
this material being in the putative dusty torus, as the dust 
sublimation radius is e stimated to be ~ 5 light months 
(Crenshaw et al. 2003|. Furthermore, our results can be 
compared with past observations of the source. For this 
we focus on the 1997 BeppoSAX observat ion of NGC 5548 
( Nicastro et al.pOOO Petrucci et al.|[2000 l and t he average 


of seven Suzaku observations conducted in 2007 (Liu et al. 


20101. In both cases, we re-fitted the archival data above 
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5 keV following the procedure described in Sect. |3.2| For 
the reflection component, we And that the spectral shape 
is always consistent with our best-fit parameters, i.e. a 
photon index of 1.9 and cut-off energy of 300 keV, while 
the reflected flux shows significant variations of about 50%. 
The primary flux (e.g., in the 2-10 keV band) shows similar 
variations, and we obtain a reflection fraction / consistent 
with that found during the 2013 campaign, i.e. / ~ 0.3, 
again calculated as the ratio between the illuminatin g and 
primary flux in the 20-40 keV range (see end of Sect. 3.21. 
Since / is consistent with being constant over a few years, 
the reflecting material should lie within a few light years 
from the primary X-ray source. In this case, the reflection 
fraction is indicative of the covering fraction of the reflector 
as seen from the X-ray source. From the definition of /, 
it follows that / = fl/(47r — fl). According to our results, 
then, the solid angle covered by the reflecting material can 
be es timated as n, consistently with previous findings 
(e.g., Nandra fc Pounds||1994 Pounds et al.||2003 |. 


Then, we have been able to constrain the primary 
continuum and the effects of absorption. The primary 
continuum is well fitted by a power law with a photon 
index varying between ^1.5 and ~1.7, with at least 
a lower limit on the presence of a high-energy cut-off 
{Ec > 50 keV), which is usually considered a manifestation 
of Comptonization and is expected to be linked to the 
coronal temperature. The approximate relation between 
the temperature kT^ of an optically thick Com ptonizing 
corona and the cut-off energy Ec is kTe = Ed'i (Petrucci 
et al. 20011. However, it should be noted that the cut-ott 
of realistic Comptonization models is very different from 
an exponential cut-off. Moreover, the estimate of the 
temperature and of the optical depth also depends on the 
geometry assumed. In one case (obs. 6) we have obtained a 
measure of 701 ^q keV for the cut-off energy, which provides 
an estimate of keV for the coronal temperature, 

with the caveats we just mentioned. Moreover, our results 
suggest a statistically significant variation of both the 
photon index and the cut-off energy. For example, in the 
case of obs. 5 we And a lower limit of 180 keV for the cut-off 
energy. We can thus estimate the coronal temperature to 
be > 60 keV for obs. 5 and < 37 keV for obs. 6, both at 90% 
confidence level. Such a variation may be due to different 
reasons, like a variation of the disk-corona geometry. The 
coronal optical depth and temperature should follow a 
univocal relationship for a given heating/cooling ratio, 
i.e. the ratio of the power dissipated in the corona to the 
intercepted soft l uminosity which cools the electrons via 
Comptoniz ation (Haardt & Maraschi 1991 Haardt et al. 
1994 19971. If the coronal optical depth changes, due for 
instance to a change in the accretion rate, the temperature 
will change accordingly, to keep the heating/cooling ratio 
constant. But the heating/cooling ratio can also change, 
responding to a variation in the disk-corona geometry, like 
a variation of the inner disk radius or a velocity variation 
of a potentially outflowing corona (Malzac et al. 20011. 
Then, the coronal optical depth and temperature will 
change in a non-trivial way to be consistent with the new 
heating/cooling ratio. Owing to these complexities, fitting 
the data with a phenomenological cut-off power law does 
not allow for a unique interpretation. However, our findings 
motivate us to apply realistic Comptonization models to 


constrain the coronal parameters. 

Concerning absorption, we needed to add two differ¬ 
ent partially covering components (only one in two cases 
out of seven). The obscurer component 1 has Nh 1 x 10^^ 
cm“^ and Cp 0.8, while the obscurer component 2 has 
Nh 9 X 10^^ cm“^ and Cp 0.3, consistently with 
the estimates of |K14[ Both show significant variability on 
a time-scale of a few days, indicating that the obscuring 
material varies along the line of sight. This relatively heavy 
obscuration was unexpected, since NGC 5548 is historically 
thought to be an archetypal unobscured Seyfert 1. The 
discussion reported in |K14| ascribes this obscuration to a 
fast and massive outflow which may be interpreted as a 
disk wind. 


Finally, we have made a first step towards a more 
physical description of the primary emission by replacing 
the simple cut-off power law with a thermal Comp¬ 
tonization model to fit the average spectrum. According 
to our analysis, the X-ray emitting corona has a mean 
temperature of 40)l^g keV and an optical depth of 2.7'dl'2i 
assuming a spherical geometry, during this campaign. We 
compared this result with past observations of NGC 5548 
by re-fitting the 1997 BeppoSAX spectrum and the 2007 
Suzaku average spectrum. We used the COMPPS model 
with spherical geometry (geom=- 4 ) like we did in Sect. 
|3.4[ The kTc — y contours are shown in fig. [m We And that 
the temperature and the optical depth are not consistent, 
at 90% confidence level, between the different observations. 
According to this test, the average values of the coronal 
temperature and optical depth show a significant long-term 
variability, namely over 10-15 years, with a clear decrease 
of the coronal temperature from about 300-400 keV in 1997 
down to 40-50 keV in 2013. The optical depth, on the other 
hand, increases from about 0.2-0.3 up to 2-3. Although 
the values quoted are specific to the assumed spherical 
geometry for the corona, the changes in temperature 
and optical depth are significant if one were to assume 
a different geometry for the corona. The corresponding 
variation of the Compton parameter y suggests a varia¬ 
tion, of uncertain origin, of the coronal heating/cooling 
ratio between these different observations. The coronal 
parameters may be variable on shorter time scales as well, 
as suggested by the detected variations of the cut-off energy. 


Comparing with the 1997 BeppoSAX observation and 
the 2007 Suzaku observations, the 2013 campaign provided 
us a broader energy range which, together with the 
high signal-to-noise ratio spectra, allows for much better 
constraints on the coronal parameters (as shown by the 
contours in Fig. [TT| . However, for the optimal study 
of ComptonizatioiqUV sp ectra may give a sign ificant 
contribution. For example, Petrucci et al. (20131 have 
been able to fit the UV and soft X-ray spectra of Mrk 
509 with a Comptonization component interpreted as a 
"warm" corona above the disk. If this scenario is also 
valid for NGC 5548, we expect a correlation between the 
UV and soft X-ray emission. This correlation can only be 
detected after the variability of the obscurer is properly 
removed. The analysis of the UV and soft X-ray emission 
by Mehdipour et al. (in prep.) indeed suggests that such 
a correlation is present, and that it may be intrinsically 
similar to that of Mrk 509. This point, and a thorough test 
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of realistic Comptonization models, will be the subject of 
a forthcoming paper. 
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